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ABSTRACT 

We report our multi-band infrared (IR) imaging of the transitional millisecond pulsar system 
J1023+0038, a rare pulsar binary known to have an accretion disk in 2000-2001. The observations 
were carried out with ground-based and space telescopes from near-IR to far-IR wavelengths. We 
detected the source in near-IR JH bands and Spitzer 3.6 and 4.5 /im mid-IR channels. Combined 
with the previously-reported optical spectrum of the source, the IR emission is found to arise from 
the companion star, with no excess emission detected in the wavelength range. Because our near-IR 
fluxes are nearly equal to those obtained by the 2MASS all-sky survey in 2000 Feb., the result indicates 
that the binary did not contain the accretion disk at the time, whose existence would have raised the 
near-IR fluxes to 2-times larger values. Our observations have thus established the short-term nature 
of the interacting phase seen in 2000-2001: the accretion disk at most existed for 2.5 yrs. The binary 
was not detected by the WISE all-sky survey carried out in 2010 at its 12 and 22 /xm bands and our 
Herschel far-IR imaging at 70 and 160 /im. Depending on the assumed properties of the dust, the 
resulting flux upper limits provide a constraint of < 3 x 10 22 -3 x 10 25 g on the mass of the dust grains 
that possibly exist as the remnant of the previously-seen accretion disk. 

Subject headings: binaries: close — stars: individual (J102347. 6+003841) — stars: low-mass — stars: 
neutron 



1. INTRODUCTION 

It was quite surprising when the binary mil- 
lisecond pulsar (MSP) J1023+0038 (hereafter 'J1023) 
was discovered in an untargeted radio pulsar survey 
(Archibald ct al. 2009J), because although th e source was 
previ ously known to have radio emission ([Bond et al.l 
2002) , it was suggested to be a cataclysmic variable (CV) 
and later has been established that it is likely a low-mass- 
X-ray binary (LMXB) with a 4.75 hour orbital period 
and a ~0.2 Mq companion ([Thorstensen fc Armstrong! 
12005ft . Detailed analyses of an optical spectrum of J1023, 
taken in 2001 February by the Sloan Digital Sky Sur- 
vey (SDSS), have shown that the source at the time 
was in a bright state and had broad, double-peaked hy- 
drogen and heliu m emission lines, a typical feature of 
an accretion disk ([Wang et al.l 12009). From 2002 May 
onward, however, the source has been seen dramati- 
cally different, having only a G-type absorption spec- 
trum ([Thorstensen fc Armst rong 2005) and thus indicat- 
ing the disappearance of the disk approximately 1.5 yrs 
later. Combining with these observational results, the 
discovery of a radio pulsar has strongly suggested that 
this is the first such binary found at the end of its tran- 
sition from a LMXB to a radio MSP. 

MSPs are believed to be "recycled" pul- 
sars, formed from evolution of X-ray binaries 
([Bhattacharva fc van den Heuvell Il991| ). In such a 
binary with a low-mass companion, a neutron star was 
spun up by accretion of matter from the companion 
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via an accretion disk. At the end of the interactive 
phase when the mass accretion rate in the disk was 
below a minimum value and thus the inner edge of 
the disk was truncated outside of the light cylinder of 
the neutron star by the neutron star's magnetic field, 
the ne utron star appeared at radio frequencies as a 
MSP ([Campana et al.l 119981 and references therein). 
Subsequently, the disk would be disrupted due to the 
radiation pressure from the fast-spinning pulsar. While 
it is generally accepted, there are few observational 
studies of the scenario due to the lack of known systems 
that arc right at the end of the transition. 

From the estimated spin-down rate of the J1023 pul- 
sar, its rotational energy loss rate (so called spin-down 
luminosity) has be en found to be ~ 4.3 x 10 34 ergs s _1 
([Deller et al.ll2012l) . The major part of the energy should 
be carried out in a magnetized, high- velocity wind, which 
is likely interactin g with the companion causing th e or- 
bital modulation ([Thorstensen fc Armstrong] I2005T) and 
eclipses of the pulsar emission ([Archibald eHH T2009) 
seen at optical and radio frequencies, respectively. The 
wind also plausibly caused the disappearan ce of the disk 
that existe d in 2000-2001 (however see iTakata et al.l 
120101 [20121 for an alternative scenario). 

It is not clear whether the companion star in J1023 
would repeatedly have an outflow to form a disk. If it 
does, as pointed out bv lThorstensen fc Armstrong! (j2005f ) 
that the Kelvin-Hclmholtz relaxation time for the com- 
panion is much lo nger than 10 years since 2001 (see also 
iDeller et al.ll2012T) . J1023 would serve as a rare example 
for detailed studies of features of the transitional pro- 
cesses and pulsar-wind-disk interaction. 

In an effort to fully study this rare system and the 
detailed evolutionary process of such binary systems at 
or right after the transition phase, we have carried out 
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infrared (IR) observations with ground-based and space 
telescopes. The accretion disk that has disappeared in 
the optical after 2002 May might have left an IR remnant. 
Recent Chandra observations of the source have found or- 
bital variability in its X-ray emission and li kely revealed 
an intrabinary shock ([Bogdanov et al.ll2011h . The shock, 
which would be produced by the interaction between the 
pulsar wind and outflow from the companion, might also 
be detectable at IR wavelengths. Our observations and 
data reduction are described in §[2j The results from the 
observations are given in § [3j and the implications of our 
results for J1023 are discussed in §|4] 

2. OBSERVATIONS 
2.1. Near-IR Observation and Data Reduction 

Our target was detecte d in near-IR JHK S b ands by 
the 2MASS all-sky survey (jSkrutskie et al.ll2006l) in 2000 
February 6, but with relatively large uncertainties (J = 
16.30 ± 0.10, H = 1 5.69 ± 0T3, and. K s = 15.9 with 
no uncertainty given; IWang et alll2009H In order to ei- 
ther confirm the 2MASS results or detect any variability, 
we re-observed J1023 in JH bands using the 2.5-meter 
Irenee du Pont telescope at Las Campanas Observatory 
in Chile on 2012 May 10. The camera was the RetroCam, 
which uses a detector of a Rockwell Hawaii-1 1024 x 1024 
HgCdTe array. The field of view was 3.4' x 3.4', with a 
pixel scale of 0.2" pixel -1 . The total on-source exposure 
time at each band was 405 s. During an exposure, the 
telescope was dithered in a 3x3 grid with offsets of ap- 
proximately 15" to obtain the measurement of the sky 
background. The observing conditions were good, hav- 
ing 0.6" seeing. 

We used the IRAF data analysis package for data reduc- 
tion. The images were dark-subtracted and flat-fielded. 
From each set of dithered images in one observation, a 
sky image was made by filtering out stars. The sky im- 
age was subtracted from the set of images, and then the 
sky-subtracted images were shifted and combined into 
one final image of the target field. 

The near-IR counterpart to J1023 was well detected 
in both J and H band images. We performed aperture 
photometry to measure brightnesses of the target and 
other in-field sources. Flux calibration was conducted 
by using the 2MASS sources detected in our images. 

2.2. Spitzer IRAC 3.6 and 4-5 fim Imaging 

We observed the target on 2010 January 15 four times 
with the Spitzer Space Telescope during its warm mis- 
sion phase (Program ID 60197). The multiple exposures 
were requested in order to average out any possible or- 
bital flux changes caused by the irradiated companion 
star by the pulsar. The imagin g instrument used was the 
Infrared Array Camera QRAC: lFlzi"o et al.ll200l . While 
the camera operated in four channels at 3.6, 4.5, 5.8, and 
8.0 /im, only the first two channels were available during 
the warm mission and used for our observations. The de- 
tectors at the short and long wavelengths are InSb and 
Si: As devices, respectively, both with 256x256 pixels and 
a plate scale of 1.2" pixel -1 . The field of view (FOV) is 
5'x5'. The frame time was 30 s with 23.6 s and 26.8 s 
effective integration time per frame, and the total inte- 
gration times in each of our four observations were 1.97 
min and 2.23 min at channel 1 and 2, respectively. 



The raw image data were processed through the IRAC 
data pipelines (version S18.13.0) at the Spitzer Sci- 
ence Center (SSC). In the Basic Calibrated Data (BCD) 
pipeline, standard imaging data reductions, such as re- 
moval of the electronic bias, dark sky subtraction, flat- 
fielding, and linearization, are performed and individual 
flux-calibrated BCD frames are produced. The details 
of the data reduction in the pipelines can be found in 
the IRAC Instrument Handbook (version 2.0.2). Using a 
MOsaicker and Point source Extractor (MOPEX) pack- 
age, provided by the SSC, we conducted further data re- 
duction and profile-fitting photometry. In the processes, 
we used the MOPEX pipelines Overlap to perform back- 
ground matching between the BCD frames, Mosaic to 
carry out removal of cosmic rays and bad pixels, and APEX 
to perform point source extraction. Also using Mosaic, 
final post-BCD (PBCD) mosaics were produced by com- 
bining the BCD frames. The target was well detected in 
our IRAC images. 

2.3. WISE Imaging 

Launched on 2009 December 14, the Wide-field In- 
frared Survey Explorer (WISE) mapped the entire sky 
at 3.4, 4.6, 12, and 22 /im (called Wl, W2, W3, and 
W4 bands, respectively) in 2010 with FWHMs of 6.1", 
6.4", 6.5", and 12 0" in the four bands, respectively 
(jWright et all IMtCl) . The WISE all-sky images and 
source catalogue were released in 2012 March. We 
checked the WISE data and found that J1023 was de- 
tected in Wl and W2 bands but not in the other 
two bands. We downloaded the images of the target 
field from the Infrared Processing and Analysis Center 
(IPAC). The dates of the observations and the depth of 
coverage of the target field were between 2010 May 21- 
26 and 11-12 pixels (corresponding to 97-106 s on-source 
integration time), respectively. 

2.4. Herschel PACS 70 and 160 fim Imaging 

We also observed J1023 on 2011 November 26 with 
the Herschel Space Observatory (Obs. ID 1342233048, 
1342233049). The imaging instrument used was the 
Photodetector Array Camera and Spectrometer (PACS; 
iPoglitsch et al.l[2010l ) . PACS performs imaging at a blue 
and a red channel simultaneously, of which the central 
wavelengths are either 70 or 100 and 160 /im, respec- 
tively. The detectors are a 64x32 and a 32x16 bolome- 
ter array for the blue and red channels, respectively, both 
with a field of view of ~ 1.75' x 3.5'. We used the mini- 
scan map mode of PACS for imaging at the channels of 
70 y, (bandwidth 25 /im) and 160 /jm (bandwidth 85 /im). 
The total telescope time for our observation was 1 hour 
with 18 min on-source time. 

The image data were downloaded from the Herschel 
Science Archive, which were processed from the Herschel 
Data Processing system at the Herschel Science Centre 
(HSC). The target was not detected in the images. We 
used the background counts at the source region and es- 
timated the flux upper limits. The values we derived 
were highly consistent with that given by the Herschel 
observation planning tool HSpot. 

3. RESULTS 

The obtained flux measurements for J1023 from 
2MASS JHK S1 our near-IR JH, and Spitzer/IRAC 3.6 
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and 4.5 /xm detections are given in Table [TJ We as- 
signed a 0.3 mag uncertainty to the 2MASS K s band 
measurement. The JH magnitudes resulting from our 
observation are consistent with the 2MASS values within 
the uncertainties. For the Spitzer observations at each 
IRAC band, the four flux measurements showed <5% 
variations, and their averaged value is listed in the table. 
In addition the WISE 3.4 and 4.6 /im flux measurements 
from the WISE catalogue are also provided. For the non- 
detections of J1023 in the WISE W3 and W4 bands and 
our Herschel images, flux upper limits were obtained and 
their 3er values are given in Table [TJ 
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Fig. 1. — IR flux measurements and upper limits of J 1023 
(triangles: near-IR JH bands; squares: 2MASS measurements; 
plus signs: Spitzer /IRAC 3.6 and 4.5 ^m mid-IR channels; di- 
amonds: WISE four mid-IR bands; X signs: Herschel/PACS 
70 and 160 )im bands). The averaged optical spectrum from 
Thorstcnscn & Armstrong (2005) is also plotted as a solid curve. 
A simple blackbody spectrum (dotted curve) well describes the 
near-lR and mid-IR flux measurements, indicating non-detection 
of excess IR emission at the wavelength range. A model spectrum 
of the accretion disk in 2001 and the combined spectrum of the 
disk and blackbody are over-plotted (dash dotted and long-dash 
dotted curves, respectively). Model spectra from two sizes of dust 
grains are shown (dashed and long-dashed curves are for 7 and 20 
fim dust grains, respectively), providing upper limits on the dust 
mass of 3 X 10 22 and 3 X 10 25 g for the first and latter size dust. 

We note that the WISE measurements of J1023 in the 
two bands, particularly W2 band, have flux values lower 
than those of the IRAC, although they are consistent 
within 3<7 uncertainties of the WISE fluxes. To compare 
the measurements of the two telescopes and check the 
reliability of the WISE measurements, we obtained flux 
measurements of a few WISE sources that were also de- 
tected in our IRAC images. The WISE measurements 
generally matches those of the IRAC down to 15 mag, 
and the differences become large for sources fainter than 
15 mag. We calculated the average differences result- 
ing from 12 sources with similar brightnesses to that of 
J1023, and they are 0.06±0.19 mag for IRAC/3.6 mi- 
nus WISE/3.4 and -0.02±0.33 mag for IRAC/4.5 minus 
WISE/4.6. As seen in Table [TJ the same differences for 
J1023 are -0.12 and -0.67 mag. Therefore the WISE 
and IRAC flux values are also consistent within the ~2cr 



systematic differences. 

3.1. Near-IR and Mid-IR detections 

Optical emission from J1023 after 2002 May is known 
to arise from the companion star, which has a G5V 
type optical spectrum and probably has a size ap- 
proximately equal to that of the Roche lobe (0.4 i? ; 
IThorstensen fc Armstrong 120051; iWang et al.ll2009D . We 
di splay the optical s pectrum (the same spectrum given 
in IWang et al.1 120091 which was averaged from 23 spec- 
tra o btained in 2003-2004 by IThorstensen fc Armstrong! 
120051 and dereddened with A v = 0.14) and IR fluxes 
in Figure 1, in whi ch the near-IR flux es are dereddened 
with A v = 0.14 (jWang et al.l 12009ft and the redden- 
ing at the mid-IR bands are ignored. It can be seen 
that the near-IR and mid-IR emission we detected also 
likely arose from the companion star. To confirm this, 
we fit the data points ( 2MASS near-IR and WISE 
band 1 and 2 data points were not included) with a 
simple blackbody model. In order to reduce the dom- 
inance of the optical spectrum, which has 1500 data 
points, we assigned a weight of 1/300 to them. We 
found that a blackbody temperature of 6000 K and a 
radius-to-distance ratio R/D of 6.6xl0~ 12 can well de- 
scribe our data (see Figure [JJ). The temperature value 
is clos e to 5700 K found bylThorstcnsc n fc Armstrong! 
(2005) from fitting the modulation light curves with an 
irradiated-star mod el. Given D = 1.37 kpc measured by 
iDeller et al.l (j2012h . R « 0.4i? Q , which is equal to the 
Roche-lobe radius. Therefore our fit results are consis- 
tent with the known properties of the binary system, and 
we conclude that no excess IR emission was found from 
our detections. 

3.2. Implication of the near-IR measurements 

It was noted by IWang et al.l (|2009l) that the 2MASS 
measurements of J1023 made on 2000 Feb. 6 were con- 
sistent with being at the tail of the G5V spectrum of the 
companion. Combining with the mid-IR measurements, 
our more accurate near-IR JH measurements, which are 
approximately equal to the 2MASS values, confirm the 
2MASS detection and the companion-star origin of the 
near-IR emission detected both in 2000 and at the cur- 
rent time. This result raises a very interesting point that 
in 2000 Feb. the accretion disk did not exist. In Fig- 
ure [TJ we over - plotte d the disk model spectrum shown 
in IWang et al.l <|2009f ) and the combined spectrum of the 
disk and blackbody. As can be seen, if the accretion disk 
had existed in 2000 Feb., the near-IR fluxes would have 
been approximately 2 times larger, J = 15.2, H = 15.0, 
a nd K„ = 14. 8 , base d on the combined model spectrum. 

iBond et al.1 ([2002) presented an optical spectrum of 
J1023, which was taken on 2000 May 6 and contains the 
same emission features as the SDSS spectrum. The exis- 
tence of the accretion disk thus likel y started after 2000 
Feb. 6 but before May 6. Previously IWang et al.l (|2009f ) 
found a required time of 0.5 yrs for building the accre- 
tion disk in 2001 by considering the estimated disk mass 
of 10 23 g and mass transfer rate of 10 16 g s . Our 
constraint on the starting time of the disk existence is 
therefore roughly consistent with their estimation. The 
life time of the interacting phase should be at most 2.5 yrs 
from 2000 Feb. to 2002 May, suggesting that for J1023 
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the interacting phase is of short-term in nature and that 
the system might repeat the short-term process in the 
future when the outflow from the companion could over- 
come the pulsar wind. Close monitoring of the binary 
system is warranted. 



4. DISCUSSION 

We have conducted multi-band imaging observations of 
J 1023 from near-IR to far-IR wavelengths using ground- 
based and space telescopes. In the near-IR JH bands 
and Spitzer mid-IR 3.6 and 4.6 /jm channels, we have de- 
tected the source but the emission has been found to arise 
from the companion star of the binary and no excess IR 
emission was detected. Presumably any hot (600-800 K) 
dust in the system that has a flux of ~0.1 mJy, which 
is approximately equal to the fluxes of the companion in 
the mid-IR bands, would have been detected. We have 
thus set such an upper limit on hot-dust emission. At 
the longer wavelengths, we did not detect emission from 
the source and the 3<r upper limits were approximately 
from 0.5 mJy at 12 /im to 6 mJy at 160 fj,m. 

Constraints on the existence of relatively cold dust 
material in the binary system can be estimated from 
our results. A remnant of the accre tion disk that was 
disrupted by either the pulsar wind (IWang et al. 12009 ) 
or 7- ray emission from the pulsar (Ta kata et al. I20K , 
2012) might stay in the system and be illuminated by 
the pulsar wind. We us ed a simple dust model given by 
iFoster fc Fischer] ((1996), in which a single temperature is 
assumed for all dust grains, to estimate infrared emission 
from possible cold dust material in J1023. The spin-down 
lu minosity of the puls ar in J1023 was recently estimated 
bv lDeller et al.l ()2012f l. Considering that 1% of this en- 
ergy output heats the dust of two sizes, 7 and 20 ^tm, 
we found that our IR upper limits constrain the masses 
of the dust to be lower than 3xl0 22 g and 3xl0 25 g, 
respectively. The corresponding dust temperatures are 
^300 and ~80 K. The model spectra are d isplayed in 
Figure [Q We note that IWang et al] (|2009D have esti- 
mated 10 23 g mass for the accretion disk seen in 2001. 
Presuming that 1% of the remnant's mass is the dust, 
the upper limits on the total mass of the remnant are 
-10 24 g for the dust of 300 K or -10 26 g for the dust 



of 80 K. The masses are substantially higher than the 
mass of the accretion disk, suggesting that the observa- 
tions were not deep enough (particularly at 8 ^m with 
no Spitzer observations) for the detection of the possibly 
existing dust material in the system. 

X-ray emission from J1023 has b een detected with a 
dominant non-thermal component ([Homer etaLl [20061: 
lArchibald etai]|2010t iBogdanov et al.ll2011lh Combined 
with the fact that orbital variability was seen in the X- 
ray emission, the origin of the emission is likely an in- 
trabinary shock, pro duced by the outflows from the two 
stars of the system (|Bogdanov et al.l [20111 : iTakata et al.l 
2012). This shock could radiate emission in a long wave- 
length range from X-ray to IR wavelengths (see, e.g., 
ISlane et all l2008f) . However, assuming a single power- 
law e mission and extending the Chandra spectrum (a = 
— 0.2: lBogdanov et al.ll201lf ) to the IR wavelengths, a flux 
at 10 /im is predicted to be 10~ 4 mJy. Therefore much 
deeper mid- or far-IR observations than ours are needed 
in order to determine possible emission properties of the 
putative intrabinary shock at the long wavelengths. 
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Facility: Irenee du Pont (RetroCam), Spitzer (IRAC), 
Herschel (PACS) 
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TABLE 1 

Flux measurements and upper limits of J1023 



Telescope /Instrument 


Observation date 


Band 


Magnitude 


Flux (mjy) 


2MASS 


2000 Feb 6 


J 


16.30±0.10 


0.481±0.045 






H 


15.69±0.13 


0.542±0.065 






K a 


15.9±0.3 a 


0.29±0.08 


Du Pont/RetroCam 


2012 May 10 


J 


16.260±0.024 


0.499±0.011 






H 


15.964±0.025 


0.421±0.019 


Spitzer/IRAC 


2010 Jan 15 


3.6 


15.870±0.034 


0.126±0.004 




4.5 


15.811±0.038 


0.085±0.003 


WISE 


2010 May 21-26 


3.4 


15.989±0.066 


0.123±0.008 






4.6 


16.49±0.34 


0.043±0.014 






12 


<12.1 


<0.4 






22 


<8.8 


<2.6 


Herschel/PACS 


2011 Nov 26 


70 




<2.8 






160 




<6.3 



Note. — For non-detections, 3c flux upper limits are given. 
0.3 mag uncertainty is assigned by us to K 3 



